Spin Dynamics in Stripe-Ordered Las/3Sri/3Ni04 
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Polarized and unpolarized neutron inelastic scattering has been used to measure the spin excita- 
tions in the spin-charge-ordered stripe phase of La5/3Sr!/ 3 Ni04. At high energies, sharp magnetic 
modes are observed characteristic of a static stripe lattice. The energy spectrum is described well 
by a linear spin wave model with intra- and inter-stripe exchange interactions between neighbouring 
Ni spins given by J = 15 ± 1.5 meV and J' = 7.5 ± 1.5 meV respectively. A pronounced broadening 
of the magnetic fluctuations in a band between 10 meV and 25 meV is suggestive of coupling to 
collective motions of the stripe domain walls. 



Over the past decade, the tendency of doped antiferro- 
magnetic oxides to exhibit symmetry-broken phases in- 
volving the ordering of both spin and charge has become 
increasingly apparent. Widespread interest was gener- 
ated by the discovery of a stripe-like, spin-charge or- 
dered phase in a non-superconducting layered cuprate [Q . 
This stripe phase, the like of which has now been found 
in many other doped transition metal oxide systems, is 
characterized by parallel lines of holes that act as charged 
domain walls separating regions of antiferromagnetically- 
ordered spins. Its discovery fuelled a debate about the 
role played by stripe correlations in the formation of the 
superconducting state in the cuprates 0, and stimu- 
lated numerous experimental investigations into stripe 
pchenomena. These have focussed principally on the 
La2- x Sr x Ni04_|-a series, which exhibits stripe ordering 
over a wide range of hole concentration |3j Q . 

While the static properties of ordered stripe phases are 
now quite well characterized, less is known about their 
dynamics. Specifically, there is a lack of information on 
how the collective motions of the holes in a stripe do- 
main wall couple to the spin dynamics of the antiferro- 
magnetic domains. One way to make progress in under- 
standing stripe dynamics is to use neutron inelastic scat- 
tering to probe the spin excitation spectrum as a function 
of wavevector and energy in simple compounds exhibit- 
ing well-correlated stripe ordering. As well as providing 
information on the microscopic interactions that govern 
the properties of stripes, such studies can address the 
question of whether stripes are essential or incidental to 
the mechanism of superconductivity. 

Here we report polarized- and unpolarized-neutron 
scattering measurements of the spin excitation spectrum 
in the ordered stripe phase of La5/ 3 Sr!/ 3 Ni04. This com- 
position has not been studied before by neutron inelas- 
tic scattering, and was chosen because the spin-charge 
order is particularly well correlated and has a very sim- 
ple superstructure commensurate with the crystal lattice 
Pj @j H> ||] — see Fig. |l](a) . Compared with earlier neu- 
tron scattering studies of stripe-ordered La2_ a: Sr x Ni04-|_ ( 5 
compounds |5|, |(| our measurements cover a wider range 
of wavevector and energy than previously explored in any 
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FIG. 1: (a) Model for the stripe order found in the Ni-0 



plane of Las/aSrj^NiCU 



Arrows denote 5 = 1 spins 



on the Ni 2+ sites, open circles represent Ni 3+ holes. The O 
sites are not shown. <f> is the angle between the spin axes 
and the stripe direction. Primitive unit cells of the N1O2 
square lattice and stripe superstructure are indicated respec- 
tively by the small square and parallelogram. J and J' are 
respectively the intra- and inter-stripe exchange interactions 
between nearest-neighbour Ni spins, (b) Diagram of recipro- 
cal space showing several Brillouin zones for the stripe-order 
superlattice shown in (a). The (h, k) indices correspond to 
the NIO2 square lattice, and the small circles are the stripe 
superlattice zone centres. The lines marked A and B indicate 
the directions of the scans used in the experiment. 



one compound, and the simplicity of the magnetic struc- 
ture of La 5 / 3 Sr 1 / 3 Ni04 allows us to perform a quanti- 
tative analysis of the data by linear spin-wave theory. 
Thus, for the first time in a stripe-ordered compound, 
we determine values for the nearest-neighbour coupling 
strengths for spins within a stripe domain and for spins 
separated by a charged domain wall. We also observe 
an anomaly in the magnetic scattering between lOmcV 
and 25meV that could, we suggest, originate from cou- 
pling between spin excitations and collective motions of 
the charged domain walls. 

Two single crystals of La 5 / 3 Sr 1 / 3 Ni04 in the form 
of rods, 7-8 mm in diameter and ~40mm in length, 
were used for the experiments. Both were grown from 
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the same batch of starting powder by the floating-zone 
method [ flO| . Their oxygen excess was determined by 
thermogravimetric analysis to be S = 0.01 ± 0.01, and 
their electrical and magnetic properties were found to be 
in very good agreement with data in the literature. 

Neutron inelastic scattering measurements were per- 
formed on the IN20 and IN22 triple-axis spectrometers 
at the Institut Laue-Langevin. The incident and final 
neutron energy was selected by Bragg reflection from an 
array of either pyrolytic graphite (PG) or Heusler alloy 
crystals, depending on whether unpolarized or polarized 
neutrons were employed. Scans were performed with a 
fixed final energy of either Ef = 14.7 meV or 34.8 meV, 
and a PG filter was present to suppress higher-order har- 
monics in the scattered beam. Two settings of the crys- 
tal were used, giving access to the (h, h, I) and (h, —2h, I) 
planes in reciprocal space. In this paper the reciprocal 
lattice is indexed with respect to a body-centred tetrago- 
nal lattice with cell parameters a — 3.8 A and c = 12.7 A. 

Long-range ordering of holes and spins in 
La5/3Sr 1/ / 3 Ni04 occurs below Tqo — 240 K and 
Tso — 200 K respectively. The stripe arrangement 
formed in the Ni-0 layers is shown in Fig.|l](a) ||, |). 
In this model the spins are collinear, and <p is the angle 
between the spin axes and the stripe direction. All the 
measurements reported here were taken at T = 14 K, 
at which temperature 4> — 53° Q and the in-plane and 
out-of-plane spin and charge order correlation lengths 
are several hundred A and 20-50 A respectively |8| . In 
inelastic measurements the c axis correlations were found 
to decay very rapidly with energy, entirely disappearing 
above 5 meV. Over the energy range considered here it 
is therefore reasonable to treat the spin dynamics as 
two-dimensional. 

Fig. 0(b) shows several two-dimensional Brillouin zones 
of the stripe order. We probed the excitation spectrum 
either by scanning the neutron scattering vector Q along 
the lines marked A and B at a fixed energy, or by scan- 
ning the energy at a fixed Q. In reality, the ordered 
stripe phase contains equal proportions of two twins with 
stripes at 90° to one another. The reciprocal space for 
the second twin is rotated by 90° with respect to that 
shown in Fig.[l](b). Experimentally, one observes a su- 
perposition of the scattering from both twins. 

Assuming the ordered ground state in Fig. [l] one ex- 
pects to observe magnetic excitations dispersing from the 
magnetic zone centres, and indeed this is confirmed by 
our measurements. To illustrate the results, Fig. || shows 
wavevector scans parallel to the (h, h, 0) and (h, —2h, 0) 
directions measured with unpolarized neutrons at a fixed 
energy of 70meV. The six peaks observed in Fig.||(a) 
correspond to spin excitations propagating in a direction 
perpendicular to the stripes away from successive zone 
centres along the A line in Fig.|l|(b). Similarly, the two 
strongest peaks in Fig. ^|(b) are spin excitations propagat- 
ing away from (1/3,-2/3) along line B, approximately 
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FIG. 2: Q scans parallel to (a) (h, h, 0) and (b) (h,-2h,0) 
at a constant energy of 70 meV, showing the magnetic scat- 
tering from La 5 /3Sr!/3Ni04 The lines are simulations of the 
scans generated by convolution of the calculated spectrometer 
resolution with the spin- wave model discussed in the text. 
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Dispersion of the magnetic excitations in 
La 5/ /3Sr 1/ / 3 Ni04 parallel to (a) (h,h,0) and (b) (h, — 2h, 0). 
The full and broken lines are the spin-wave model dispersion 
for the two twins, calculated with parameters J = 15meV, 
J' = 7.5 meV, and K c = 0.07 meV. 



parallel to the stripes. Under the experimental condi- 
tions used we were not able to resolve the separate peaks 
associated with each zone centre when the energy was 
below 37meV. 

In Figs. H we plot the dispersion of the spin excitations 
along the (h, h, 0) and (h, —2h, 0) directions as deduced 
from a series of (mainly) constant-energy scans like those 
of Figs.||. To arrive at the data points shown in Figs.|| 
we corrected the observed peak positions for small shifts 
(10-20%) caused by the spectrometer resolution as it 
intersects the curved dispersion surface. We determined 
the corrections by convolution of the resolution function 
with the model dispersion surface to be described later. 

In the energy range 10-30 meV very strong scatter- 
ing from phonons makes it difficult to identify unam- 
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FIG. 4: (a) Spin-flip (SF) and non-spin-flip (NSF) scat- 
tering as a function of energy at a constant wavevector of 
(4/3,4/3,0). The sharp peaks in the NSF channel are due 
to phonons. The SF background is estimated from data col- 
lected at (1.2,1.2,0). (b) SF scattering intensity plotted as 
a function of wavevector for energies of 10 meV, 15 meV and 
26 meV. The upper two scans have been shifted vertically by 
150 and 300 counts respectively. The lines are fits to Gaussian 
peak profiles. A monitor of 5 x 10 corresponds to a counting 
time of 5-8 mins, depending on energy. 



biguously the magnetic scattering with unpolarized neu- 
trons. Therefore, to explore this energy range we em- 
ployed polarization analysis, keeping the neutron polar- 
ization direction P parallel to Q by adjusting currents 
in a Hclmholtz coil mounted around the sample position. 
In this configuration the neutron spins are flipped when 
scattered by magnetic excitations, but remain unchanged 
in non-magnetic (e.g. phonon) scattering processes. 

Fig.[|(a) shows the neutron spin-flip (SF) and non- 
spin-flip (NSF) scattering observed as a function of en- 
ergy with Q fixed at the (4/3, 4/3, 0) magnetic zone cen- 
tre. The separation of the magnetic and non-magnetic 
scattering by the polarization analysis is exemplified by 
the sharp phonon peaks centred on 16.5 meV and 22 meV 
in the NSF channel that are not present in the SF chan- 
nel. However, the most remarkable feature in this scan is 
the broad valley centred near 15meV in the SF scatter- 
ing. By comparison with the SF measurements made at 
Q = (1.2, 1.2,0), where the magnetic scattering has de- 
creased to background level, we see that the amplitude 
of the magnetic signal above background is nearly a fac- 
tor 3 smaller at E = 15meV than at either E = lOmeV 
or E = 26meV. Similar scans measured at (1/3,1/3, 1), 
I = 4.5, 5.5, show the same valley feature. A secondary 
feature in Fig.^(a) is the fall in magnetic intensity with 
energy observed below ^7meV. 

One possible origin for the features just described could 
be spin anisotropy gaps in the excitation spectrum. Be- 
low an anisotropy gap a component of the spin fluctu- 
ations freezes out, causing a reduction in scattering in- 
tensity. We checked this possibility by comparing the 
magnetic scattering intensities for neutron polarizations 



P || Q and P _L Q. The data at lOmeV, 15meV 
and 26meV were all consistent with the assumption of 
isotropic transverse spin fluctuations about the direction 
of the ordered moment. At 3meV, however, the signal 
is predominantly due to in-plane transverse spin fluctua- 
tions. Therefore, we conclude that the drop in intensity 
below ~7meV is caused by an out-of-plane anisotropy 
gap, but that the main valley feature cannot be explained 
by spin anisotropy. 

Fig. |](b) shows the SF scattering intensity measured in 
a series of constant-energy scans in the neighbourhood of 
the valley feature. These data show that, in addition 
to a reduction in amplitude, there is also an increase in 
width. The 15meV peak is found to be ~50% broader 
in wavevector than the lOmeV and 26meV peaks. This 
broadening suggests that the valley feature is associ- 
ated with the decay of antifcrromagnetic spin excitation 
modes through hybridisation with another type of exci- 
tation. One possibility is a coupling to optic phonons 
intrinsic to the material, but we do not believe this to be 
the case because (a) the valley region is ^10meV wide, 
much larger than the width of individual phonon exci- 
tations, and (b) we do not observe any structure in the 
SF scattering channel due to the excitation of phonons 
through a magnetoelastic interaction. Instead, we sug- 
gest that a probable explanation for the valley feature is 
the interaction between spin excitations and motions of 
the stripe domain walls with characteristic frequencies in 
the 3-5 THz range. 

Having established the overall features of the spin ex- 
citation spectrum we now describe the spin-wave model 
used to extract the microscopic exchange parameters. 
The model is based on the spin Hamiltonian 



H 



z\2 



(1) 



where the first two summations are over pairs of nearest- 
neighbouring Ni spins, the first sum having both spins 
within the same stripe domain and the second having 
the two spins in adjacent domains separated by a line 
of holes. J and J' are the corresponding exchange pa- 
rameters — see Fig. 0(a). The third term describes the 
out-of-plane anisotropy. The in-plane anisotropy is much 
smaller than the out-of-plane anisotropy, and so is ne- 
glected. 

We calculated the energy dispersion and scattering 
cross-section using the linear spin-wave approximation, 
and after convolving the model spectrum with the spec- 
trometer resolution we compared the results with the 
measured scans. The parameter values giving the best 
agreement with the totality of data are J = 15 ± 1 .5 meV, 
,/' = 7.5 ± 1.5 meV and K c = 0.07 ± 0.01 meV. For 
reference, the J and K c parameters obtained by Naka- 



jima et al |11 from a fit to the spin- wave dispersion of 
La2Ni04 were 15.5 meV and 0.52 meV respectively. We 
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find, therefore, that the spin anisotropy is strongly re- 
duced in the doped compound, but the strength of the 
nearest-neighbour exchange interaction in the antiferro- 
magnetic region is essentially unaffected by hole-doping. 

The lines drawn on Figs.Ba) and (b) illustrate scan 
simulations, and are seen to match the main features of 
the data quite well. The inclusion of more exchange 
parameters might achieve even better agreement with 
experiment, but is not justified given the extent of the 
present data set. We did, however, consider an alterna- 
tive model in which J' couples pairs of spins displaced 
by (a, a) across a domain wall instead of (2a, 0). How- 
ever, the deviations from the first model were slight and 
mainly confined to energies close to the zone boundary 
energy where the data is sparse. Therefore, while there 
may be systematic errors due to the linear approximation 
and the neglect of more distant couplings, we expect that 
the present analysis yields a good estimate of the relative 
strengths of the intra- and inter-stripe spin couplings. 

The results presented here provide several new insights 
into the dynamics of stripe phases. At high energies 
(> 30 meV) the spin excitations propagate as if the un- 
derlying charge-ordered lattice were static. A similar 
observation has been made by Bourges et al who 
also commented on a lack of significant anisotropy in the 
spin-wave velocity in their study of Lao.i.69Sro.3iNi04. 
The advantage of the La5/ 3 Sr 1 / 3 Ni04 compound studied 
here is its simple stripe superstructure that permits us 
to quantify the anisotropy in the spin- wave dispersion in 
terms of a microscopic model. Thus, while our data are 
consistent with those of Bourges et al, our analysis has 
determined the ratio J/ J' to be close to 2 and uncov- 
ered a dramatic reduction in single-ion spin anisotropy 
relative to La2Ni04. 

Concerning the anomalous broadening observed in the 
energy range 10-25 meV, we earlier discussed several 
pieces of evidence that lead us to suggest that this broad- 
ening arises from a coupling between spin excitations and 
the collective motions of the stripe domain walls. Further 
to these, we note that the stripe ordering temperature 
(~240K) in this compound translates to ^20mcV. In 
other words, the thermal energy needed to destroy long- 
range order of the stripes is comparable to the energy 
where we observe the anomalous broadening in the spin 
excitations. Zaanen et al fl2]| discussed the dynamics 
of charged domain wall motion in an antiferromagnetic 
background on general theoretical grounds, and argued 
that domain wall fluctuations will have their strongest in- 
fluence on the spin excitation spectrum at low energies. 
In one reading this is consistent with our results, but 
it is interesting that the anomalous broadening appears 
to be restricted to a band of energies below which the 
spin excitations recover their sharp profile. A possible 
explanation is that there exists a commensurability gap 
of ^lOmeV for domain wall motions due to the pinning 
of the charges to the lattice. If this were the case, then 



one might expect the valley feature to be absent from 
the spin excitation spectrum of compounds whose stripe 
period is incommensurate with the lattice. 

Finally, we draw attention to recent calculations of the 
imaginary parts of the charge and spin dynamical sus- 
ceptibilities of an ordered stripe system described by the 
Hubbard model (p| Q. In addition to transverse spin 
waves these calculations also predict longitudinal modes 
arising from meandering and compressive movements of 
the domain walls. We did not observe sharp phason-likc 
modes of the type predicted, but their weak intensity rel- 
ative to the spin wave scattering and the broadening in 
the valley region may have precluded their observation. 
Similar calculations for the special case of 1/3 doping, in- 
corporating the experimentally-determined exchange pa- 
rameters, would be valuable for interpreting the features 
we have measured, and would thus provide a direct test 
of this theoretical description of stripe dynamics. 
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